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Abstract 

We investigate the particle trajectories in an irrotational shallow water flow 
over a flat bed as periodic waves propagate on the water's free surface. Within 
the linear water wave theory, we show that there are no closed orbits for the 
water particles beneath the irrotational shallow water waves. Depending on the 
strength of underlying uniform current, we obtain that some particle trajectories 
are undulating path to the right or to the left, some are looping curves with a 
drift to the right and others are parabolic curves or curves which have only one 
loop. 

1 Introduction 

The motion of water particles under the waves which advance across the water is 
a very old problem. The classical description of these particle paths is obtained 
within the framework of linear water wave theory. After the linearization of the 
governing equations for water waves, the ordinary differential equations system 
which describes the particle motion turns out to be again nonlinear and explicit 
solutions are not available. In the first approximation of this nonlinear system, 
one obtained that all water particles trace closed, circular or elliptic, orbits (see, 
for example, [13], [16], [20], [21], [22], [23]), a conclusion apparently supported 
by photographs with long exposure ([32], [22], [23]). Consequently there is no 
net transfer of material particles due to the passage of the wave, at least, at this 
order of approximation. 

While in these approximations of the nonlinear system all particle paths appear 
to be closed, in [10] it is shown, using phase-plane considerations, that in linear 
periodic gravity water waves no particles trajectory is actually closed, unless 
the free surface is flat. Each particle trajectory involves over a period a back- 
ward/forward movement, and the path is an elliptical arc with a forward drift; 
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on the flat bed the particle path degenerates to a backward/forward motion. 
Similar results hold for the particle trajectories in deep-water, that is, the tra- 
jectories are not closed existing a forward drift over a period, which decreases 
with greater depth (see [5]). These conclusions are in agreement with Stokes' 
observation [23]: "There is one result of a second approximation which may 
possible importance. It appears that the forward motion of the particles is not 
altogether compensated by their backward motion; so that, in addition to their 
motion of oscillation, the particles have a progressive motion in the direction 
of the propagation of the waves. In the case in which the depth of the fluid is 
very great, this progressive motion decreases rapidly as the depth of the particle 
considered increases." 

For shallow water waves, the standard results are that the orbits described by 
water particles beneath waves are elongated ellipses with the longer axis paral- 
lel to the flat bottom, and at the bottom the orbits are straight lines (see, for 
example, [T5]V 

Similar conclusions hold for the governing equations without linearization. 
Analyzing a free boundary problem for harmonic functions in a planar domain, 
in [4] it is shown that there are no closed orbits for Stokes waves of small or large 
amplitude propagating at the surface of water over a flat bed; for an extension of 
the investigation in [3] to deep-water Stokes waves see [TB] ■ Within a period each 
particle experiences a backward/forward motion with a slight forward drift. In a 
very recent preprint [3] , the results in [3] are recovered by a simpler approach and 
there are also described all possible particle trajectories beneath a Stokes wave. 
The particle trajectories change considerably according to whether the Stokes 
waves enter a still region of water or whether they interact with a favorable or 
adverse uniform current. Some particle trajectories are closed orbits, some are 
undulating paths and most are looping orbits that drift either to the right or to 
the left, depending on the underlying current. 

Analyzing a free boundary problem for harmonic functions in an infinite planar 
domain, in [B] it is shown that under a solitary wave, each particle is transported 
in the wave direction but slower than the wave speed. As the solitary wave 
propagates, all particles located ahead of the wave crest are lifted while those 
behind have a downward motion. 

Notice that there are only a few explicit solutions to the nonlinear governing 
equations: Gerstner's wave (see [14] and the discussion in 2])), the edge wave 
solution related to it (see [3]), and the capillary waves in water of infinite or finite 
depth (see [T2], [IS]). These solutions are peculiar and their special features 
(a specific vorticity for Gerstner's wave and its edge wave correspondent, and 
complete neglect of gravity in the capillary case) are not deemed relevant to sea 
waves. 

The present paper is concerned with the particle trajectories in an irrota- 
tional shallow water flow over a flat bed as a periodic wave propagates on the 
water's free surface. It is natural to start this investigation for shallow water 
waves by simplifying the governing equations via linearization. In Section 2 we 
recall the governing equations for water waves. In Section 3 we present their 
nondimensionalisation and scaling. The linearized problem in the irrotational 
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shallow water regime is written in Section 4. We also obtain the general solution 
of this problem. The next section is devoted to the description of all the possible 
particle trajectories beneath a linear periodic irrotational shallow water wave. 
We see that these particle trajectories are not closed. Depending on the strength 
of underlying uniform current, denoted by the constant en, we obtain that: for 
c > 2 the particle trajectories are undulating path to the right, for cq < — 1 the 
particle trajectories are undulating path to the left, for — 1 < en < the particle 
trajectories are looping curves with a drift to the right and for < Co < 2 the 
particle trajectories are parabolic curves or curves which have only one loop. 



2 The governing equations for gravity water waves 

We consider a two-dimensional inviscid incompressible fluid in a constant grav- 
itational field. For gravity water waves these are physically reasonable assump- 
tions (see [TB] and [H]). Thus, the motion of water is given by Euler's equations 

u t +uu x + vu z = -- p Vx 
v t + uv x + vv z = - jPz - g 

Here (x, z) are the space coordinates, (u(x, z, t),v{x, z, t)) is the velocity field of 
the water, p{x, z, t) denotes the pressure, g is the constant gravitational acceler- 
ation in the negative z direction and p is the constant density. The assumption 
of incompressibility implies the equation of mass conservation 

u x + v z = (2) 

Let ho > be the undisturbed depth of the fluid and let z = ho + i]{x, t) repre- 
sent the free upper surface of the fluid (see Figure 1). 



A 




z=0 Bottom 



Figure 1 . A periodic water wave propagating over a flat bed. 



The boundary conditions at the free surface are constant pressure 

p = pa on z = h + i](x, t), (3) 

Po being the constant atmospheric pressure, and the continuity of fluid velocity 
and surface velocity 



v = rj t + u-q x on z = ho + rj{x, t) 



(4) 
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On the flat bottom z — 0, only one condition is required for an inviscid fluid, 
that is, 

» = oiiz = (5) 

Summing up, the exact water-wave problem is given by the system {I}-©. In 
respect of the well-posedness for the initial- value problem for (HJ-© there has 
been significant recent progress, see and the references therein. 

A key quantity in fluid dynamics is the curl of the velocity field, called 
vorticity. For two-dimensional flows with the velocity field (u(x, z, £), v(x, z, t)), 
we denote the scalar vorticity of the flow by 

u(x,z) = u z - v x (6) 

Vorticity is adequate for the specification of a flow: a flow which is uniform 
with depth is described by a zero vorticity (irrotational case) , constant non-zero 
vorticity corresponds to a linear shear flow and non-constant vorticity indicates 
highly sheared flows. 

The full Euler equations (HJ-© are often too complicated to analyze directly. 
One can pursue for example a mathematical study of their periodic steady 
solutions in the irrotational case (see Q], [25]) or a study of their periodic steady 
solutions in the case of non-zero vorticity (see [7J, [5])- But in order to reach 
detailed information about qualitative features of water waves, it is useful to 
derive approximate models which are more amenable to an in-depth analysis. 



3 Nondimensionalisation and scaling 

In order to develop a systematic approximation procedure, we need to char- 
acterize the water-wave problem ([l])-© in terms of the sizes of various fun- 
damental parameters. These parameters arc introduced by defining a set of 
non-dimensional variables. 

First we introduce the appropriate length scales: the undisturbed depth of wa- 
ter ho, as the vertical scale and a typical wavelength A (see Figure 1), as the 
horizontal scale. In order to define a time scale we require a suitable velocity 
scale. An appropriate choice for the scale of the horizontal component of the 
velocity is y/gho- Then, the corresponding time scale is ^ h and the scale for 

the vertical component of the velocity is ho . The surface wave itself leads 
to the introduction of a typical amplitude of the wave a (see Figure 1). For 
more details see [16]. Thus, we define the set of non-dimensional variables 

x i-> Ax, z i-> h z, rj^ai], 1 1-> -j=j^t, 
u i ^ \fgha~u, v i-» ho ^ a £° v 

where, to avoid new notations, we have used the same symbols for the non- 
dimensional variables x, z, rj, t, u, v, on the right-hand side. The partial 
derivatives will be replaced by 
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Let us now define the non-dimensional pressure. If the water would be station- 
ary, that is, u = v = 0, from the equations ([1]) and ([3]) with rj — 0, we get for 
a non-dimensionalised z, the hydrostatic pressure po + pgho(l — z). Thus, the 
non-dimensional pressure is defined by 

P t-^ p + pgh (l - z) + pgh p (9) 

therefore 

Px p—Px, Pz -pg + pgpz (10) 

Taking into account ©, ©, © and ([TO)), the water-wave problem ([J)-© 
writes in non-dimensional variables, as 

u t + uu x + vu z = -p x 
5 2 (v t + uv x + vv z ) = -p z 

u x + v z = (11) 
v = e(r/ t + urj x ) and p = eq on z = 1 + tr](x, t) 
v = on z = 



where we have introduced the amplitude parameter e = and the shallowness 
parameter 6 = fy". In view of ([§]). the vorticity equation © writes in non- 
dimensional variables as 

u z — 5 2 v x + ^®^° lj(x, z) (12) 
9 

For zero vorticity flows (irrotational flows) this equation writes as 

u z = 6 2 v x (13) 

After the nondimensionalisation of the system ©-© let us now proceed 
with the scaling transformation. First we observe that, on z = 1 + erj, both 
v and p are proportional to e. This is consistent with the fact that as e — > 
we must have v — > and p — ?> 0, and it leads to the following scaling of the 
non-dimensional variables 

p^?ep, (u, v) h-t e(u, v) (14) 

where we avoided again the introduction of a new notation. The problem (|11[) 
becomes 

u t + e(uu x + vu z ) = —p x 
S 2 [v t + e(uv x + vv z )\ = -p z 

u x + v z = (15) 
v — r/t + e.w] x and p — i] on z = 1 + £T]( x i t) 
v = on z = 
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and the equation (fT2"|) keeps the same form. 

The system which describes our problem in the irrotational case is given by 



u t + e(uu x + vu z ) = -p x 
8 2 [v t + e(uv x + vv z )} = -p z 
u x + v z = 



u z = S 2 v r 



(16) 



rjt + e.w] x and p — r\ on z = 1 + a](x, t) 
v = onz = 



4 The linearized problem 

The two important parameters e and S that arise in water-waves theories, are 
used to define various approximations of the governing equations and the bound- 
ary conditions. The scaled version (fT6"|) of the equations for our problem, allows 
immediately the identification of the linearized problem, by letting e — > 0, for 
arbitrary d. The linearized problem in the shallow water regime is obtain by 
letting further 5 — > 0. Thus, in the irrotational case, we get the following linear 
systems 

u t + Px = 
Pz -0 

u x + v z = , 
u z =0 {l) 
v = rjt and p = r\ on z = 1 
v = on z = 

From the second equation in (TTTl) we get in the both cases that p does not 
depend on z. Because p = rj(x, t) on z = 1, we have 

p = rj(x,t) for any < z < 1 (18) 

Therefore, using the first equation and the fourth equation in (|17[) . we obtain 
in the irrotational case 

t 

rj x (x, s)ds + J-(x) (19) 

o 

where T is an arbitrary function such that 

F(x)=u(x,0) (20) 



Differentiating (|19p with respect to x and using the third equation in (|17l) we 
get, after an integration against z, 

v = -zu x = z (^J T] xx (x, s)ds - T' (x) J (21) 
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In view of the fifth equation in (fT?)) we get after a differentiation with respect 
to t, that rj has to satisfy the equation 

Vtt -r)xx=0 (22) 

The general solution of this equation is r)(x,t) — f(x — t) + g(x + t), where / 
and g are differentiable functions. It is convenient first to restrict ourselves to 
waves which propagate in only one direction, thus, we choose 

V (x,t)=f(x-t) (23) 

From (|2Tj) , (|23p and the condition v = r\t on z = 1, we obtain 

T(x) = f(x) + c (24) 

where cq is constant. 

Therefore, in the irrotational case, taking into account (TT51) . (fig)) . (|2ip . (|23|) and 
(j2"4")) . the solution of the linear system (IT71) is given by 

V {x,t)=f{x-t) 

p(x,t) = f(x-t) , 25 > 

u(x, z, t) = f{x - 1) + c 

f (x, z, £) = —zf'{x — t) = —zu x 

5 Particles trajectories in the irrotational case 

Let (x(t), z(t)) be the path of a particle in the fluid domain, with location 
(x(0), z(0)) at time t = 0. The motion of the particle is described by the 
differential system 



dx - u(x,z,t) 
§ = v(x,z,t) 



itZinVK ( 26 ) 



with the initial data (a;(0),z(0)) := (a; ,zo)- 
Making the Ansatz 

f{x -t)= cos(27r(a; - t)) (27) 
from (|25|) . the differential system (|26|) becomes 



=cos(2vr(a;-t))+co 

(28) 

^ = 27rzsin(27r(a; - t)) 

Notice that the constant Co is the average of the horizontal fluid velocity over 
any horizontal segment of length 1, that is, 

1 f x+1 

c = - I u(s,z,t)ds, (29) 
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representing therefore the strength of the underlying uniform current. Thus, 
Co = will correspond to a region of still water with no underlying current, 
Co > will characterize a favorable uniform current and Co < will characterize 
an adverse uniform current. 

The right-hand side of the differential system (|2"5)l is smooth and bounded, 
therefore, the unique solution of the Cauchy problem with initial data (xq, zq) 
is defined globally in time. 

To study the exact solution of the system (|2"5|) it is more convenient to re- 
write it in the following moving frame 

X = 2ir(x-t), Z = z (30) 

This transformation yields 

^ = 2ttcos(X) +27r(c - 1) 

(31) 

4f = 27rZsin(X) 
Let us now investigate the differential system (|3"TT) . 

5.1 The case c = 

The horizontal component of the velocity u in (|25l) . with f(x — t) given by (|27|) . 
has in the moving frame (|30j) . the following expression 

u{X,Z,t) =cos(X) + c (32) 

Thus, the case Co = is obtained for 

u(X,Z,t)dX = (33) 

This is the Stokes condition for irrotational flows, that is, the horizontal velocity 
has a vanishing mean over a period. 

In the considered case, we write the first equation of the system (|3"Tj) into the 
form 

dX , N 

cos(a ) - 1 

We use the following substitution (see [17], 1.76, page 308) 

2v v 2 — 1 2 

sin(X) = ^^, cospO = ^— , dX = -^—dy (35) 

y* + 1 y + 1 y + 1 

In the new variable, integrates at 

y = 2nt + k (36) 
k being an integration constant. Hence, 

X(t) = 2arccot (2nt + k) (37) 
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Taking into account (|55|) . we obtain 



sva(X(t)) 



2(2irt + k) 



1 + (2nt + k) 2 
Therefore, the second equation in (f5T|) yields 



(38) 



Z(t) = Z(0)exp / 27rsin(X (s))ds = Z(0)exp In 



1 + (2nt + k) 2 



l + k 2 



From (13T)|) . (|3T|) and (IMj) . we obtain that the solution of the system 
the initial data (x , Zq), has the following expression 

x(t) = t + iarccot (2nt + k) 
z{t) = T ^[l + {2^t + kf] 

From the initial conditions, we get k := cot (tvxq). 

The derivatives of x(t) and z(t) with respect to t, have the expressions 



(39) 
with 



(40) 



, ( ,S _ (27Tt + k) 2 -l 

Z '(t) = ^(2nt + k) 



Therefore 



x'(t) > 
z'{t) > 



\2wt + k\ > 1 
{2-irt + k) > 



the flat bottom being at z = 0, we have zq > 0. 

Thus, for tin the intervals (-co, =^), {=^ k , -£;), (-^r, ) and ( ' °°)' 
the derivatives cc'(i), z'(i), have the following signs 



t : 





-l-k 




1-fe 




2tt 




2tt 


x'{t) > 


x'(t) < 


< 


x'(t) > 


a/(t) < 


z'(t) < 


z'[t) > 


z'(i) > 



(41) 



The limits of x(t), z(t) and ^& for t — ^ — oo and i — ?• oo are 



lim x(i) 



lim x(t) — oo 



(42) 



Thus, taking into account (1411) and (J42J) , we sketch below the graph of the 
parametric curve (|40[) 



lim z(t) — oo, lim z{t) = oo 

t— > — oc i— >oo 
r *(*) r 

lim — — - = — oo, lim — — = oo 

i^-oc x{t) t^o° X(t) 
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Figure 2. Particle trajectory in the case c Q =0 

Thus, we get: 

Theorem 1: In the case of no underlying current, the particle trajectories be- 
neath the irrotational shallow water waves are curves which have only one loop 
like in Figure 2. 

5.2 The case c (c — 2) > 

In this case, we write the first equation of the system (f3"TT) into the form 

/ 1 ^ = 27rf (43) 

J cos(X)+c -l V ' 

We use the same substitution (j35j) . In the new variable, (l43l) becomes 

dy 



= 2nt (44) 



which integrates at 



, ^— arctan (. °° y] = 2nt + k (45) 
co V c o ~ 2 VV c o - 2 y ) V 1 

k being an integration constant. Further, we obtain 

y = Cotan(a(t)), (46) 

with 

*o := (47) 
V c 

a(t) :=-^-{2Ttt + k) (48) 
Hence, returning to the variable X, we get 

X(t) = 2arccot \<L tan (a(t)) 1 (49) 
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Taking into account (|55|) . (jlS)) , we obtain 

2£ tan (a(t)) 



sin(X(t)) 



1 



£q tan (<*(*)) 



(50) 



The second equation in (|3ip yields 



= Z(0)exp ( y 2tt sin(X (s))ds 



(51) 



From ([50]). gSJ) and J5TJ, we obtain that the solution of the system ([55]l, with 
the initial data (xq,Zq), zq > 0, has the following expression 



x(t) = t+ ^arccot <£q tan (a(t)) 



(52) 







2 


1+ 


Co tan(a(s))) 





The derivatives of x(t) and with respect to t, have the expressions 

_ (C,j D.-i.,--^,/,,)...-, 1 _ 2 ( S in 2 (a( t ))-^— 



' T ^ cos 2 (a(t)) [l+Cg tan 3 (ct(t))j 



z '( t ) = ZQ 4^£ tan( a (t)) 



1 + 



C tan(o(t)) 



' ex P I Jo — : 



(-co) 



Co tan(a(s)) 



cos 2 (a(t))[l+Cg tan 2 (a(t))J 



Co tan(a(s))) 



yds 



Let us now study the signs of the derivatives in (|53|) . We are in the case 
cq(cq — 2) > 0, that is, cq € (— oo, 0) U (2, oo). 



(a) If c < -1, then > L Therefore, sin 2 (a(i)) 



< 0. Thus, 



we obtain that x'(t) < 0, for all t. 
The sign of z'(t) will depend on a(t). For a(t) in intervals of the form a(t) € 
(— f + i7r, In), I £ Z, we get z'(i) < 0, and for a(t) € (/7r, Z-zr + §) , Z G Z, we get 
> 0. 

We sketch below the particle trajectory in this case: 



Figure 3. Particle trajectory in the case c < - 1 
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(b) If 1 < c < 0, then < 1. Thus, 



for a(t) < - arcsin (^f^^j + In, we get x'(t) > 0, z'(t) < 0, 



for a(t) e - arcsin J ^£ + ?tt, In , we get ic'(i) < 0, z'{t) < 0, 



for a(t) G Ztt, arcsin V 2^ + /tt we get cc'(t) < 0, z'(t) > 0, 



for a(t) > arcsin I y ° 2 J + In, we get x'(t) > 0, > 0, 
where / € Z. We sketch below the particle trajectory in this case: 




Figure 4. Particle trajectory in the case - 1<. c Q < 

(c) If c > 2, then ^_£i > 1. Therefore, sin 2 (a(t)) - < q. Thus, we 

obtain that x'(t) > 0, for all t. 

The sign of z'(t) will depend on a(t). z'(t) < for a(t) <G (— f + /7r, In), I € Z, 
and > for a(t) G (/tt, ivr + f ) , Z G Z. 

We sketch below the particle trajectory in this case: 



Figure 5. Particle trajectory in the case c > 2 



Therefore, we proved: 

Theorem 2: In the case that the underlying uniform current is moving in the 
same direction as an irrotational shallow water wave and the strength of the 
current is bigger than 2, then the particles trajectories beneath the wave are 
undulating paths to the right (see Figure 5). 

In the case that the underlying uniform current is moving in the opposite 
direction as an irrotational shallow water wave and the strength of the current is 
smaller than -1, then the particles trajectories beneath the wave are undulating 
paths to the left (see Figure 3). If the strength of the adverse current is bigger 
than -1, then the particle trajectories are loops with positive drift (see Figure 
4). 
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5.3 The case c e (0, 2] 

In this case, we write the first equation of the system (|5T1) into the form 

J cospO+co-1 
We use the same substitution (1351). In the new variable, (l54l) becomes 



dy 



c o J y 



2-c 



2vrf 



(54) 



(55) 



which integrates at 



C() 



c V 2 - c 



■In 



2 -en 
co 



2-co 
co 



2-Kt + k 



(56) 



k being an integration constant. Further, we obtain 

exp(2j0(t)) + 1 



y = &o 



where 



exp(2/3(t)) - 1 
exp(2/3(f)) - 1 ; 
exp(2/3(f)) + 1 

£o := 



if \y\ > #0 , 

if \y\ <*<>, 



2 -co 



CO 

/3(f) :=^(27rf + fc) 
Hence, returning to the variable X, we get 



(57) 

(58) 
(59) 



X(t) = 2arccot 



JSo coth (/3(f)) or X(t) = 2arccot j?o tanh (/3(f)) (60) 



if | cot (y) I > -So, respectively, | cot (=§-) I < -So 
Taking into account (|55|) . ([57| , we obtain 



sin(X(f)) = 



2£q coth (/3(f)) 



coth (/3(f)) 



sin(X(f)) = 



2^o tanh (/3(f)) 



1 



J?o tanh (/3(f)) 



(61) 

Thus, the solution of the system (|28|). with the initial data (:ro, Zq), z n > 0, has 
in this case the following expressions 



x(t) =t + iarccot coth (/3(f)) 



z(f) - z exp /„* ^coth(^ M ) , 

\ 1+ Ro coth(/3(s))) / 



(62) 
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or 



x(t) = t + -arccot tanh (/?(£)) 
z(t) = z Q exp ( /„* 

\ 1+ £0 tanh(/3(s))) 



rd,S 



We derive x{t) and z{t) from ([62]) with respect to i and we get 



(63) 



= 1 + 



ah 2 (/5(t)) 



2 -co 



1+.S 2 C oth 2 (0(t)) 



= ZQ 4^ocoth(ff(*)) exp Jt 4^ coth(/^)) 



1 + 



8.0 coth(/3(t))) 



1 + 



.So coth(/3(s))) 



(64) 



Because we are in the case Co € (0, 2], we have 2 — cq > 0. Thus, the derivative 
a/(i) > for all t. The sign of z'(t) depends on the sign of f3(t), that is, for 
0(t) < we have z'(t) < and for /3(i) > we have z'(t) > 0. Then, the 
particle trajectory in this case is like in Figure 6 (a). 
For the second alternative ([63]) . we get 



x'{t) 



(^ + l)sinh 2 (/3(t))+c Q -l _ 
Dsh 2 (/3(t))[l+^ 2 tanh 2 (,3(t))J _ 



2 <BixHor[p^t))-- 



cosh 2 (,3(t))[l+.^ tanh 2 (/3(t))J 



_ ZQ 4^q tanh(/3(s)) gxp / jt 4irfj tanh(/3(s)) 



1 + 



£0 tanh(/3(s))) 



ji tanh(/3(s))) 



(65) 

(a) If 1 < c < 2, then we get x'(t) > 0, for all t. 

The sign of z'(t) will depend on (3(t). For fj{t) < 0, we get z'(t) < 0, and for 
/3(f) > 0, we get z'(t) > 0. 

The particle trajectory in this case is like in Figure 6 (a). 



(b) If < c < 1, then 



for smh(/3(fj) < - J ^2., we get x'(t) > 0, z'{t) < 0, 



for sinh(/3(i)) € - 



go— eg 



, we get x'(t) < 0, z'(t) < 0, 



for sinh(/3(i)) S 0, 



^ we get ir'(t) < 0, z'(t) > 0, 



for sinh(/3(t)) > y /£ V a > we § et > °> > °> 

Thus, the particle trajectory in this case is sketched in Figure 6 (b). 



5 Particles trajectories in the irrotational case 



15 




Figure 6. Particle trajectory in the case < Cg < 2 



We thus have: 

Theorem 3: In the case that the underlying uniform current is moving in the 
same direction as an irrotational shallow water wave and the strength of the 
current is smaller than 2, then the particles trajectories beneath the wave are 
parabolic curves or curves which have only one loop like in Figure 6. 
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